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Abstract 

> ■ 

00 ■ A semi-phenomenological model is proposed to study dynamics and 

[ — ' steady-states of magnetic fluxes and currents in mesoscopic rings and 

. cylinders at non-zero temperature. The model is based on a Langevin 

equation for a flux subject to zero-mean thermal equilibrium Nyquist 
noise. Quenched randomness, which mimics disorder, is included via the 
fluctuating parameter method. In the noiseless case, the stability thresh- 
old (critical temperature) exists below which selfsustaining currents can 
run even if the external flux is switched off. It is shown that selfsustaining 

S and persistent currents survive in presence of Nyquist noise and quenched 
disorder but the stability threshold can be shifted by noise. 

^ ! PACS number(s): 05.40.-a, 73.23.Ra, 02.50.Ey 

o 

1 Introduction 



Quantum phenomena manifested at the mesoscopic level have attracted much 
5_j ■ experimental and theoretical attention. Phase coherence and persistent currents 

can be mentioned as examples. Persistent currents were predicted as early as 
1938 JjJ and have been observed experimentally only since 1990 ||. In the paper 
we analyze the steady states of magnetic fluxes and currents in a mesoscopic 
system subject to dissipation and fluctuations ||. The system consists of a set 
of concentric one dimensional rings stacked along a certain axis. It is known 
that a small metallic ring threaded by a magnetic flux displays a persistent cur- 
rent, signifying quantum coherence of electrons (called coherent electrons) and 
it is a direct manifestation of phase-coherent electron motion and the Aharonov- 
Bohm effect at the mesoscopic level. Moreover, it has been shown Q that in 
such a system selfsustaining currents run even if the external flux is switched 
off. At temperature T = 0, the system is in the ground state and only co- 
herent electrons exist. Then the persistent current flows without dissipation. 
At temperature T > the amplitude of the persistent current run by coherent 
electrons decreases and some electrons become "normal" (i.e. non-coherent). 
The motion of normal electrons is random and their flow is dissipative. Under 
some conditions, coherent conduction and normal conduction coexist, resulting 
in dissipation of a total current. It was confirmed experimentally in Q that 
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mesoscopic rings connected to a current source presented an Ohmic resistance 
which was not zero. 

We introduce a semi-phenomenological model describing the flux and current 
in such a system, which takes into account both dissipation and fluctuations. 
The approach is based on the notion of the two fluid model [|| of normal and 
coherent electrons formulated as a Langevin equation with a noise term. In the 
system at temperature T > there are various sources of noise and fluctuations. 
There are so-called universal conductance fluctuations j?J that arise from the 
random quantum interference between many electron paths which contribute to 
the conductance in the diffusive regime. These fluctuations decay algebraically 
with temperature and can be neglected at higher temperatures fl. Inelastic 
transitions in the ring cause another kind of fluctuations. However, they do 
not destroy persistent currents but decrease their amplitude There is also 
a part of the current noise which is called shot noise [|| , the spectral density of 
which is proportional to mean current. This noise can be reduced by increasing 
the size of rings Thermal motion of charge carriers in any conductor is a 
source of random fluctuations of current which is called Nyquist noise || . This 
thermal equilibrium noise is universal and exists in any conductor, irrespective 
of the type of conduction and is analogous to noise driving the motion of Brow- 
nian particles. Moreover, this noise increases with temperature. Therefore at 
relatively high temperatures and relatively large rings, universal conductance 
fluctuations and shot noise can be neglected, and only Nyquist noise can play 
an important role. This is the case considered in the paper. We assume that co- 
herent and normal electrons coexist || and that normal electrons are subjected 
to Nyquist noise. This noise generates the flux fluctuations which indirectly 
influence persistent currents run by coherent electrons. Our main goal is to 
answer the question whether persistent and selfsustaining currents survive in 
the presence of dissipation and fluctuations. In classical physics dissipation and 
fluctuations may often be described by introducing random forces into evolution 
equations. In quantum physics, the inclusion of fluctuations requires more care 
because quantum systems are described by Hamilton operators and it ensures 
conservation of energy. In the paper we combine a classical Langevin equation 
with noise term and with terms of a quantum origin. Our model is minimal in 
the sense that in two limiting cases it reduces to the well-established models of 
the quantum persistent current of coherent electrons and the classical Nyquist 
current of normal electrons. The approach used could be justified in a more 
elegant way applying the methods of the thermo-field dynamics . 

The paper is structured as follows. In Sec. II, we describe the model based on 
the concept of two fluids and introduce the Langevin equation for the magnetic 
flux, which include Nyquist fluctuations. Its detailed analysis is carried out 
in Sec. III. In this section, influence of quenched disorder on the system is 
analyzed as well. Conclusions are made in Sec. IV. Finally, in the Appendix, 
we present the Nyquist relation from which the intensity of thermal fluctuations 
is determined. 

2 Two fluid model of electronic transport 

Persistent currents in conducting structures with topology of the ring are an 
example of a quantum-interference phenomenon in which the phase coherence 
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of the electron wave function affects equilibrium and transport properties of 
mesoscopic systems. At zero temperature, T = 0, the current in a quasi 
one-dimensional mesoscopic ring threaded by the magnetic field of the Bohm- 
Aharonov type is paramagnetic |Ti"| 

Ieven(<t>) = h V" A„ Sm ( ] (1) 

in the case of an even number of coherent electrons and diamagnetic ]TT[ ] 

Iodd{4>) = Ieven{4> + <t>o/2) ( 2 ) 

in the case of an odd number of coherent electrons. Here, <f> is the magnetic 
flux, 0o := h/e is the flux quantum, the amplitudes 

2 

A n = — cos(nkplx), (3) 
irn 

the unit current I — heN e /(2l^,in), kp is the Fermi momentum, l x is the 
circumference of the ring and N e is the number of coherent electrons in the 
ring. 

Now, let us consider a collection of rings, individual current channels, which 
form a cylinder. There are N z channels in direction of the cylinder axis and 
N r in the direction of the cylinder radius. We assume that the thickness of 
the cylinder wall is small if compared with the radius. Because of the mutual 
inductance between rings, the current in one ring induces flux in other rings. 
In turn, the flux induces a current, and so on. We will analyze the effect of the 
mutual inductance among the rings. We assume that the rings are not contacted. 
So, there is no tunneling of electrons among the channels and the charge carriers 
moving in the different rings are independent. It has been shown in |l^ ] that the 
effective interaction between the ring currents, when taken in the selfconsistent 
mean field approximation, results in the magnetic flux <j) — LI to t felt by all 
electrons, where L is the cylinder inductance and Itot is the total current in a 
cylinder. For a cylinder of the radius r and the height l z the inductance is fl3|| 

2 

irr 

L = fJ- -j-, (4) 

where /io is the permeability of the free space. At temperature T > 0, the 
current I co h(4>,T) of the coherent electrons in a set of N — N r x N z current 
channels forming the cylinder is either Q 

Icoh{<t>, T) = I even {<t>, T) = NI Y A n (T) sin ( (5) 

for an even number of coherent electrons in each single channel or 

Icoh{4>,T) = I dd{4>,T) = I e , e „(0 + 0o/2,T) (6) 
for an odd number of coherent electrons. The amplitude 

a /m\ 4T exp(-nT/T*) , 

A - {T) - ^ i-exp^r/V 03 ^- (7) 
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The characteristic temperature T* is given by the relation fc^T* = Af/2tt 2 , 
where fcs is the Boltzmann constant and Ap is the energy gap at the Fermi 
surface. For temperatures T < T* the coherent current flows in such a cylinder 
without dissipation but its amplitude ^ is reduced On the other hand, 

at temperature T > 0, normal electrons occur and their flow is dissipative. The 
motion of normal electrons is random, like the motion of electrons in a normal 
conductor and it generates random currents. 

The current coming from the normal electrons can be induced by e.g. the 
change of the magnetic flux <f>. According to the Lenz's rule and the Ohm's law 
one gets Q 

RInortt) = (8) 

where R is the effective resistance of the system ||. 

At temperature T > 0, coherent and normal electrons can coexist resulting 
in dissipation of the total current Itot which, in the absence of fluctuations, is a 
sum of the persistent and normal currents, 

Itat{<t>,T) = Iooh(<P,T) +/„or(0). (9) 

The magnetic flux is related to the total current via the expression 

<p = (j> ex t + LI tot {<P,T), (10) 

i.e. it is a sum of the external flux <fi ex t and the flux coming from the currents. 
Combining (||), (g) and (|l0|) yields the equation 

~ = ™fa-&*t) + W6T)- (ii) 

It is known that the current-flux characteristics for the coherent electrons is 
extraordinary sensitive to a change of parity of the coherent carriers number 
JlT| . In order to take into account the possible difference of parity in the 
rings we consider the current of coherent electrons as the average ^cokC^'-^) = 
pleven{4>i T) + (1 — p)I dd{4 ) i T) , where p £ [0, 1] is the probability of the even 
number of coherent electrons in a given channel. 

The term describing current fluctuations should also be included in the evo- 
lution equation ( pi] ) . Within our model, the only source of fluctuations is the 
equilibrium Nyquist normal current noise generated by the resistance R. The 
correlation function of this source of fluctuations is assumed to be given by the 
Nyquist relation. It follows from the Appendix that Eq. ( pi] ) should be modified 
to the following form 



~ = - + ™« + sj ^ T(t), (12) 

where T(t) is Gaussian white noise (see the Appendix). This equation takes the 
form of a classical Langevin equation and is our basic evolution equation. 

Let us introduce dimensionless variables. In the Langevin equation (|l^), the 
basic quantity is the magnetic flux <f> — <p(t) . The natural unit of the flux is the 
flux quantum 0o = h/e. Accordingly, the flux is scaled as x — 4>/4>q. To identify 
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the characteristic time To, let us consider a particular case of (|12|), namely, when 
the persistent current and the external flux are zero. Then 

R 



-£ = --4> + \/2RkgT r(t). (13) 
at L 

From this equation it follows that the mean value 

(<Kt)) = (0(0)) expH/T ), (14) 

where 

r = L/R (15) 

is the relaxation time of the averaged normal current. Therefore time is scaled 
as t = t/ro- In this case, Eq. ( [l2| ) can be transformed into its dimensionless 
form 

x = -V'(x) + V2D f(t), (16) 

where the dot denotes a derivative with respect to the rescaled time t and the 
prime denotes a derivative with respect to x. The generalized potential 

V(x) = V(x, A, i , p, T) = ^x 2 - Xx - i F(x,p, T), (17) 

where A = 4>ext/4>o is the rescaled external flux. The prefactor iq = NLIq/^q is a 
coupling constant characterizing the interaction between ring currents (it is the 
rescaled amplitude of the flux created by the current - it leads to selfsustaining 
currents). The function 



F(x) = F(x,p,T) = J f(x,p,T)dx (18) 

characterizes the coherent electrons and 

f(x,p,T) = pf even (x,T) + (1 - p)fodd(x,T), (19) 



where 



and 



sm(2nirx) (20) 



fodd(x,T) = f even (x + 1/2.T). (21) 



The dimensionless intensity D of rescaled Gaussian white noise T(t) = ^/tq Tfat) 
is a ratio of thermal energy to the elementary energy stored up in the inductance, 

D = h B T/e , e := ^. (22) 

The rescaled Gaussian white noise T(t) has the same moments as in (|37|). Let 
us notice that the resistance R does not occur explicitly in the rescaled equation 

H). 
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Let us evaluate the order of magnitude of the parameters appearing in our 
equations. The rescaled coupling constant 

= /i e 2 iV N e 
87rm e l z 

We assume that the cylinder has the radius r = 3-10 4 A and the height l z = 100 A. 
It consists of a set of N ~ 50 current channels in a wall of width much smaller 
than the radius. If the number of electrons in each channel is N e ~ 2 • 10 5 then 
io ~ 1. The energy gap at the Fermi surface = Ti 2 N e /(2m e r 2 ) gives the 
rescaled noise amplitude 

k B T* fi e 2 N e 



2eo 167r 3 m e l z 



(24) 



For the above values of parameters the diffusion coefficient D ~ 0.001T/T*. 
Below, unless stated otherwise, the parameter are fixed so that io = 1 and 
D = 0.001T/T*. 



3 Analysis 



In this section the properties of system described by Eq. (16) are analyzed. We 
consider in details three special cases. In the deterministic - noiseless - case, we 
neglect the influence of Nyquist noise. It is a justified approximation for very 
small intensity of noise. In the second case, we include Nyquist noise, i.e. we 
analyze the full version of (|l6|). In the third case, we extend the model assuming 
that the coupling parameter io is a random variable. It can mimic uncorrelated, 
quenched disorder of the rings. 

As follows from ([Tof), the current Itot is linearly related to the magnetic flux 
<j) (or the rescaled flux x). In a consequence, the properties and behavior of 
the current are identical to the properties and behavior of the magnetic flux. 
Therefore, below we discuss mostly the dependence of the magnetic flux on the 
parameters of the model. 

3.1 Noiseless case 

First, let us consider the deterministic case of the Langevin stochastic equation 
(ft6|) formally neglecting the Nyquist noise term f(t), i.e., 



x = -V'(x). (25) 

The stationary solutions x s of (Pq) , for which x s = 0, correspond to extrema of 
the generalized potential (jlTn, 

V'(x s )=x s -\-iof(x s ,p,T) = 0. (26) 



The solutions x s of the gradient differential equation (25) are stable provided 
they correspond to a minimum of the generalized potential ( p"7| ) and they are un- 



stable in the case of a maximum 16 . In the following we investigate properties 
of solutions x s with respect to four independent parameters: the temperature T, 
the coupling constant io which characterizes the mean-field interaction between 
rings, the probability p of the occurrence of the channel with an even number 
of coherent electrons and the external flux A. 
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Figure 1: Dimcnsionless generalized potential V(x) given by Eq. (17) is shown 
as a function of the dimensionless magnetic flux x for several values of the scaled 
temperature T/T*. The scaled amplitude io = 1 and scaled external magnetic 
flux A = 0. 

3.1.1 T and io-dependence 

The dependence of the potential (|l7| ) on the temperature for X = 0, io = 1 and 
the probability p = 1/2 is shown in Fig. 1. 

In high temperatures, only one stable solution, corresponding to zero station- 
ary flux x s — and zero current, exists. If temperature decreases, a bifurcation 
occurs - the potential becomes bistable and two non-zero symmetric minima 
appear at x s = ±x m . They correspond to two stable stationary solutions. 
Physically, it means that below some critical temperature T c the spontaneous 
flux fl7| appears and non-zero stationary current flows in the system. This 
critical temperature T c is defined by the condition that V"(x s = 0) = 0. The 
corresponding diagram is shown in Fig. 2. 

The phenomenon is analogous to the continuous phase transition in macro- 
scopic systems, and appears here as a result of the interaction of ring currents. 
The central maximum x s = xm — 0, corresponds to the unstable stationary 
solution of (p5|). 

More generally, one can notice that the stationary solutions occur where the 
linear part x — A of ( p6[ ) is equal to its periodic part iof(x,p,T). In the limit 
of io — * (very small, or no interaction of ring currents), regardless T, the 
only stationary solution of ( p6[ ) is the external flux x s = A. For intermediate io 
(typical interaction of mesoscopic rings) two stable non-zero stationary states 
can exist below T c and this number of solutions is preserved in the limit T — > 0. 
As one can infer from (p^)-(pl|), decreasing temperature enhances the periodic 
part of (|26|) but only to a maximal value defined by T — > 0. Further enhancement 
of the periodic part is possible only by increasing the coupling constant io. As 
a result of that, the critical temperature T c increases with io- Therefore, if io is 
sufficiently large (very strong interaction of rings), even more stationary states 
can occur. The number of stationary states below T c and for p — 1/2 can, in 
general, be equal to 4k — 1, {k = 1,2...) but only 2k of them of stable states. 
Lowering the temperature below T c results then in a cascade of bifurcations. The 
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Figure 2: Bifurcation of the stable stationary magnetic flux x s determined by 
Eq. (26) is depicted vs temperature T/T* for several values of i and fixed 
external flux A = 0. 



first bifurcation takes place at T — T c . With further lowering the temperature 
at T — T Cl < T c two additional pairs of stationary solutions appear and so on. 
There is one metastable and one unstable solution in every pair. The metastable 
solutions correspond to the so called flux trapped in the cylinder. Notice that in 
the limit T — > and typical io > there are always spontaneous flux solutions 
whereas the flux trapped solutions can be obtained only for sufficiently large io- 

3.1.2 The p-dependence 

In the following part of this section the temperature is set below T c . If the prob- 
ability p = 1 we have an even number of coherent electrons and paramagnetic 
current in each channel. The potential possesses two minima corresponding to 
spontaneous fluxes. Decreasing p, the probability 1 — p of finding odd chan- 
nels with diamagnetic currents increases and the spontaneous flux solutions x s 
decrease to coalesce finally into a single absolutely stable solution at x s = 0. 
The ratio p at which the coalescence occurs decreases with decreasing tempera- 
ture. Now, for sufficiently large io, five stationary states exist. Note that apart 
from the stable fluxless solution x s = there are two metastable solutions at 
5 < \x s \ < 1 and two unstable solutions. The metastable solutions correspond 
to the flux trapped in the cylinder. In realistic devices they are hardly accessible 
due to the value of the necessary parameters. Further the discussion is limited 
to the case when p = 1/2. 

3.1.3 The A-dependence 

The dependence of the stationary solutions x s on the external flux A is shown 
in Fig. 3. 

There are three different types of the generalized potential. First is a sym- 
metric double well potential which appears for A = k/2 with integral k. The 
stable solutions then always around the external flux < \x s — A| < ^. 

For the values of A close but not equal to k/2 the solutions remain in that range 
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Figure 3: Generalized potential (17) for several values of the external flux A, 
fixed i = 1 and T/T* = 0.2. 



but the double-well potential becomes asymmetric - one of the stable solutions 
becomes metastable. For the values of external flux far from half integer values 
k/2 one obtains the potential with only a single stable solution. All the men- 
tioned types of potentials are accessible for < A < 1/2 indicating a kind of the 
'structural periodicity' with respect to the external flux. 

An interesting feature of the x — A characteristic is the occurrence of the 
hysteresis loop (Fig. 4). 

With increasing A, at its certain value, the system undergoes discontinuous 
jump of x. Decreasing then the value of A, the opposite jump of a; occurs at lower 
A producing a hysteresis loop. It is a hallmark of the first order phase transition. 
The transition can occur only below the critical temperature T c . Due to the 
'structural periodicity' the hysteresis loop is repeated with the period A = 1/2 
what results in the formation of a family of loops. 



3.2 Model with Nyquist noise 

Noise and fluctuations are ubiquitous in real systems and idealization of the 
noiseless systems is sometimes not justified. In the following, we will focus 
on the system (jig) subjected to Nyquist noise. From the mathematical point 



of view, the Langevin equation (16) defines a Markov diffusion process. Its 



probability density p(x, t) obeys the Fokker-Planck equation in the form 

J=rt*. *) = •j&v'WpI?* *) + D f ) ( 27 ) 

with the natural boundary condition ]im.i x i^. 0o p(x,t) = 0. The stationary solu- 
tion p s (x) is asymptotically stable Jl^ | and takes the form 

Ps {x) = N e- v ^/° (28) 

with a normalization constant 

iV 1 = / e - v ^/D dx . (29 ) 
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Figure 4: The stationary magnetic flux x s determined by Eq.(26) vs external 
flux A and two distinct temperatures. For T/T* = 0.5 the part of the graph 
with the negative slope corresponds to unstable x s . The dependence of x s on A 
is in this case hysteretic. The amplitude io = 1. 



Let us first consider the case of absence of the external flux, A = 0. If in the 
noiseless case the system possesses only one stationary solution x s = 0, the 
probability density ( p8| ) has maximum at x — and the mean value of the flux 
(x) = 0. If in the noiseless case the system possesses three stationary states, 
the probability density ( p8| ) has three extremal points: two symmetric maxima 
which correspond to the spontaneous fluxes and one minimum at x = which 
corresponds to the unstable stationary state (see Fig. 1). Because the potential 
is reflection-symmetric, V(x) = V(—x), the mean value of any odd function 
of the flux is zero. In particular, the mean value of the flux (x) — and the 
mean value of the current is zero as well. From this point of view, properties 
of stationary states are trivial and non-zero fluxes and currents are impossible. 
However, in some situations the statistical moments are not good characteristics 
of the system because much information is lost when an integration is performed 
calculating the statistical moments pof . The relevant quantity is a stationary 
probability distribution which contains much more information about the sys- 
tem. Is any reasonable method to determine the critical value of temperature 
T c in this case? One possibility is to define the phase transition in the following 
way |2^, ^l| : the phase transition point is a value of the relevant parameter 7 of 
the system at which the profile of the stationary distribution function changes 
drastically (e.g. if a number of maxima of the distribution function changes) 
or if a certain most probable point Xq begins to change to an unstable state. 
In some cases, it is indeed a good 'order parameter' of the system. For ex- 
ample, from the measurements of the laser experiment (see e.g. p2|), one can 
obtain the stationary probability distribution of the laser intensity and one can 
observe a phase transition according to the above definition. In the case con- 
sidered here, for sufficiently low temperatures, thermal fluctuations are small 
and one expects the experimental results to be accumulated around the most 
probable values of the stationary probability distribution. It follows from ( p8| ) 
that the most probable values of the flux corresponds exactly to the stationary 
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Figure 5: Plot of the characteristic escape time r e and decay time Td as a 
function of temperature T/T* for iq = 1, A = and the strength of Nyquist 
noise D = 0.001T/T*. 

states ( p6| ) of the system (p5|). In this sense, the properties of the system are 
the same as discussed in the previous subsection. We want to emphasize that 
it is correct for low temperatures because then the residence time in a stable 
state is long. For higher temperature T, thermal fluctuations become larger. In 
turn, fluctuations of the magnetic flux around the most probable value become 
larger and larger and the residence time in a stable state becomes shorter. One 
can guess that the spontaneous current should vanish at temperature To which 
is lower than the critical temperature T c in the noiseless case. This is because 
of influence of Nyquist fluctuations. The argumentation is the following. If the 
potential is multistable then one can introduce characteristic time scales of the 
system. The first characteristic time Td — l/V"(x s ) describes decay within the 
attractor x s — ii m of the potential V(x). The second characteristic time is 
the escape time r e from the well around ±x m . This time is related to the mean 
first passage time from the minimum of the potential to the maximum. If these 
time scales are well separated, i.e. if r e >> Td then the description based on the 
most probable value seems to be correct. Otherwise, this description fails and 
we should characterize the system by averaged values of relevant variables. The 
exact formula for the escape time r e is known (2^| but is not reproduced here. 
We have calculated r e for the transition from the left minimum of the potential 
( P~T| ) assuming that the left boundary at — oo is reflecting and the right boundary 
at the maximum xm of the potential is absorbing. 

In Fig. 5 we show the dependence of two characteristic times r e and Td upon 
the rescaled temperature T/T* . The escape time r e is a monotonically decreas- 
ing function of temperature while the decay time Td monotonically increases 
with temperature. In the noiseless case, for values of parameters chosen in Fig. 
5 and from Eq.(]T7|) we estimated the critical temperature T c w 0.88T*. One 
can observe that roughly for temperatures T < 0.8T*, the characteristic time 
Td is more than one order of magnitude less than r e . Both time scales are well 
separated and selfsustaining currents are long-living states. In this sense, they 
are not destroyed by Nyquist noise. 
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Figure 6: Probability density (28) for selected values of temperature, fixed 
i = 1, A = and D = 0.001T/T*. 



Temperature T enters into the Langevin equation ( |l6| ) in two-fold: it con- 
tributes to the intensity D = ttBT/eo of Nyquist noise and to the effective 
interaction between the rings via the potential ([l7]) describing coherent elec- 
trons current. The influence of temperature on the properties of equilibrium 
distribution is plotted in Fig. 6. 

As stated before the position of the most probably values correspond to 
stable stationary states. The statistically averaged current vanishes since (x) — 
0. The stationary flux variance or mean-squared deviation a = (x 2 ) — (x) 2 = 
(x 2 ) is a non-monotonic function of temperature (Fig. 7): For T — the 
variance a — x 2 , where x s is a stationary solution of (p5[). As the temperature 
increases, a diminishes attaining a minimal value at some temperature 7\. The 
temperature T\ seems to be always larger than T c what is confirmed in the 
numerical studies. 

A further increase of temperature leads to an increase of the variance. In the 
high temperature limit, the dependence is linear as for the Gaussian distribution. 
Indeed, below the critical temperature, the distribution ( |28| ) possessing two 
peaks is clearly non-Gaussian. However, for higher temperatures the probability 
density is one-peaked. For this case, the kurtosis 



measures the relative flatness of the distribution fl28|) to the Gaussian distribu- 
tion. The kurtosis is negative and it means that the distribution ( f28[) is flat. 
It approaches zero in the high temperature limit and then the distribution (pq) 
approaches the Gaussian distribution. 

The behavior of the second moment (x 2 ) has a simple explanation in terms 
of the average energy stored in the magnetic field, i.e 



where eo is given in (|2^). For low temperatures, fluctuations are small and the 
main contribution to the energy comes from the deterministic part (f) 2 /2L. He- 



Kurt = 




- 1 



(30) 



(E) = {4> 2 )/2L = e {x 2 ) 1 



(31) 
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Figure 7: Averaged magnetic energy (E)/cq given by Eq.(31) vs scaled tem- 
perature for two values of io, fixed A = and D = 0.001T/T*. 



cause the magnetic flux <p decreases as temperature increases (cf. Fig. 2), hence 
(E) decreases as well. On the other hand, for high temperature the stationary 
probability density approaches the Gaussian distribution and in consequence the 
main contribution to the average magnetic energy comes from thermal energy, 
(E) oc kT which obviously increases when T grows. The competition between 
these two mechanisms leads to the minimal value of (E) for a certain value of 
temperature T\. 

The influence of the external field on the properties of the stationary density 
( p8| ) may be deduced from Fig. 3. Finally, let us consider the limit of a very 
weak coupling between the ring currents corresponding to a very small value 
of io. The selfsustaining stable solutions are non accessible. The solutions of 
Eq.(^5|) correspond then to the persistent currents driven by the external field. 
The stationary density forms a family of one peak curves with the most probable 
values given by A. We conclude that even in the weak coupling limit the presence 
of Nyquist noise does not destroy the persistent currents. 



3.3 Model with quenched disorder 

There are several sources of disorder which can be described as quenched fluc- 
tuations. The first one is caused by impurities distributed randomly in the 
cylinder. The presence of impurities leads to modification of the coherent cur- 
rent amplitude [^5|. Another source of disorder is of a geometrical origin. The 
radius of rings being not exactly of the same value may change randomly from 
one ring to the other. Also the number N r of current channels in the direction of 
the cylinder radius can be, in general, different in every point of vertical axis of 
the cylinder. These sources of disorder can be included into our model assuming 
that deviations of the coherent current amplitude from the average are small 
and 

io = Jo + e£, (32) 

where jo is an averaged value of the amplitude io, e characterizes intensity of 
quenched fluctuations, £ is a zero- mean random variable of values in the interval 
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Figure 8: Probability density (33) for the system with (e = 0.25) and without 
(e = 0) quenched disorder for T/T* = 0.7, j = 1, A = and D = 0.001T/T* . 

[—1,1] and of the probability density P$(z). 

The stationary probability density p s {x) of the flux is now expressed as 

p s (x) = J p(x\z)P$(z) dz, (33) 
where the conditional probability distribution 

P (X\Z) = C (z) e -lV(*)-ezF( X )]/D ( g 4) 

has the same form as (^q) with the replacement = jo + ez in the potential 
(F~^1) - Now, the normalization constant depends on z, namely, 

/oo 
e - [V ( X )-ezF(a>)]/D ^ ^ 
-oo 

In Fig. 8 we show the distribution ( |33"1) for p = 1/2, A = and the uniform 
probability distribution P^(z). The observed shift of the most probable value of 
the flux is rather small and can be interpreted as an increase of the magnitude of 
the selfsustaining current. The magnitude of peaks does not change significantly 
as well. One observes that both the depth of the minimum and the height of the 
maxima decrease. One concludes that the quenched disorder does not drastically 
change the properties of the flux with only one exception. The probability 
density in or near the critical temperature T c may change qualitatively from 
one peak to two peaks. In the critical region the quenched disorder lowers the 
critical temperature T c . 

4 Summary 

Persistent and selfsustaining currents are beautiful manifestation of quantum 
coherence in mesoscopic systems. The natural question is how do they behave 
in the presence of dissipation and fluctuations. Assuming the two fluid model 
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for mesoscopic system we have investigated the influence of Nyquist noise and 
quenched disorder on these currents. Our discussion is limited to stationary 
states of the magnetic flux and current although the proposed model of the flux 
dynamics can be, in principle, applied to study time dependent problems. 

The deterministic system, obtained as a formal noiseless limit, exhibits crit- 
ical behavior such as bifurcations (phase transitions) and hysteresis. There is a 
critical temperature T c below which non-zero selfsustaining currents flow. For 
the noisy system, we have investigated the properties of the flux probability. 
We have compared the statistical properties of the stationary probability den- 
sity with the deterministic case. Maxima of the stationary probability density 
correspond exactly to the deterministic stable solutions while the local minima 
to the unstable solutions. In the presence of randomness, there is a problem 
of defining the critical temperature. We have presented three points of view 
which are based on the mean value of the flux (no phase transitions), on the 
maxima of the probability density (the same critical properties as in the deter- 
ministic case) and on the comparison of two characteristic times (the critical 
temperature lower than in the deterministic case). 

The interesting finding is that the mean magnetic energy as a function of 
the temperature is not monotonic - for low temperatures it is decreasing while 
for high temperature it is an increasing function. The minimum is not at the 
critical temperature T c but at temperature T\ > T c . Moreover, we have included 
quenched disorder and concluded that the system is stable with respect to such 
a perturbation, although the most probable values of the stationary probability 
density increase. It might be interpreted as an enhancement of the current under 
the influence of the quenched disorder. Therefore it serves as one more example 
of the constructive role of fluctuations. The possibility of the noise-induced 
current in mesoscopic rings has been considered in J24|, ^5|. The presented 
approach can be applied not only to analysis of currents in mesoscopic cylinders 
but also in superconducting rings. There is also another class of mesoscopic 
systems where our model can be applied - carbon nanotubes. The possibility of 
persistent currents in carbon nanotubes has been investigated in |2^| . 

The general conclusion is that persistent and selfsustaining currents survive 
in the presence of the above mentioned fluctuations. However, if the intensity 
of Nyquist noise or quenched disorder is sufficiently strong they lead to the 
lowering of the temperature below which the system is in the ordered state. For 
smaller noise intensity the influence of fluctuations on coherent currents is very 
small and it is very favorable for the experimental observations of persistent and 
selfsustaining currents. 



5 Appendix 

For the paper to be self-contained, we remind one of the form of the fluctuation- 
dissipation theorem and the Nyquist relation exploited in our basic Eq. (|l2|). 
The Brownian motion of a particle of mass m in a fluid of temperature T is 
described by a Langevin equation |^7j . According to the fluctuation-dissipation 
theorem p7[ |, its form for the velocity v — v(t) reads 

mv + jv = y/2-/k B TT(t), (36) 
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where a dot denotes a derivative with respect to time, 7 is the friction coefficient, 
ks is the Boltzmann constant and T(t) is the zero-mean and Dirac ^-correlated 
Gaussian stochastic process (white noise), 

(T(i))=0, (T(t)T(s))=S(t-s). (37) 

Mutatis mutandis, the Langevin equation for the current / = I(t) in the RL 
circuit takes the form p3| 



LI + RI = y/2Rk B T r(t) . (38) 
It is one of the form of the Nyquist relation. In the case when 

<t> = LI (39) 

it can be rewritten as 



1 dd> 1 2k B T 

T(t) (40) 



Rdt L V R 
which justifies the prefactor of the noise term in Eq. (|12| ) 
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